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An evaluation of the background due to atmospheric muons in the ANTARES high energy neutrino telescope 
is presented. Two different codes for atmospheric shower simulation have been used. Results from compar- 
isons between these codes at sea level and detector level are presented. The first results on the capability of 
ANTARES to reject this class of background are given. 

1. Introduction 

The detection of neutrinos from galactic/extragalactic sources has enormous potential for the understanding of 
astrophysical processes and for an improved exploration of the far Universe. Neutrino astronomy is assuming 
an important role in experimental physics research. The principle of detection in an underwater neutrino tele- 
scope is based on the emission of Cherenkov light by muons emerging from neutrino interactions inside and 
around the detector 

The ANTARES Collaboration Q] is building an undersea neutrino telescope off the French coast at a depth of 
2500 m. The telescope will consist of 900 photomultipliers (PMTs) arranged in 12 strings. The arrival times 
of Cherenkov photons on the PMTs will be registered with an accuracy of about 1 ns, which will result in an 
angular resolution of a few tenths of a degree for the reconstruction of a high energy muon track. In principle, 
selecting upward going muons would discriminate the muon tracks induced by neutrinos (both atmospheric and 
astrophysical). In practice, however, hits due to Cherenkov light emitted by the downward going atmospheric 
muons can be erroneously reconstructed as upgoing tracks. 

A complete Monte Carlo simulation has been performed to study the response of the detector to atmospheric 
muons and to estimate the possible contamination from misreconstructed tracks, starting from the production 
of cosmic ray induced atmospheric showers, propagating the surviving muons through water to the detector, 
simulating the Cherenkov light production and the hits on the PMTs and, finally, reconstructing muon tracks. 
Atmospheric shower production and propagation through the atmosphere has been performed using two dif- 
ferent codes and the results have been compared. 

2. Atmospheric muon simulation. 

The first step of the simulation is the production of the atmospheric showers induced by primary cosmic rays. 
Two different programs have been used: HEMAS |2| and CORSIKA v.6014 |3 1. 

The HEMAS code was developed for the simulation of the atmospheric muon flux in the MACRO |4| un- 
derground detector, at the Laboratori Nazionali del Gran Sasso. It takes into account all the main physical 
processes occurring in the atmosphere: it computes the first interaction point of primary cosmic rays on the ba- 
sis of the input cross sections, propagates the electromagnetic and hadronic components of the showers, takes 
into account the deflection of charged particles due to the geomagnetic field and the Earth's curvature, allowing 
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the calculations at large zenith angle. Each hadronic interaction in the atmosphere is handled with the hadronic 
interaction code DPMJET v.II |5 1. Some approximations in the particle transport along the atmosphere restrict 
the application of the code to particle energies greater than 500 GeV. 

The CORSIKA package is the most widely used tool to model extensive air showers (HAS). It performs a 
detailed simulation of the evolution in the atmosphere of an EAS initiated by various cosmic ray particles up 
to energies of some 10^° eV. Single particles are tracked, accounting for energy losses, deflection due to mul- 
tiple scattering and the Earth's magnetic field, decays of unstable particles and electromagnetic and hadronic 
interactions. The CORSIKA code provides the option to choose one of several hadronic interaction models. In 
this work the QGSJET |6 1 model has been chosen to treat high energy hadronic interactions. 
To represent the primary energy spectrum, a simplified version of the composition model described in 1 7 1 was 
used. All the elements heavier than iron have been neglected and the remaining ones have been grouped in five 
mass groups: i) protons; ii) helium; Hi) CNO; iv) Mg; v) Fe and heavy nuclei. For each nucleus, the production 
has been subdivided into several energy intervals ranging from 1 to 2 x 10^ TeV and 2 zenith intervals: i) 
0° < t? < 60°; ii) 60° < < 85°. The total number of simulated showers is about 7 x 10^ for HEMAS 
and 2 x 10® for CORSIKA. The corresponding detector livetime varies depending on the energy and zenith 
intervals chosen, going from some hours in the lower energy range to years in the higher. 




12 3 4 5 6 7 -1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 



log,o(EyGeV) cos 

Figure 1. Atmospheric muon flux arriving at the detector surface. Left: Muon energy spectrum. Right: muon angular 
distribution. The continuous line corresponds to the HEMAS simulation, the dotted line to the CORSIKA simulation 

Muon propagation through water is performed with the MUSIC (MUon Simulation Code) program fS^ . MUSIC 
is a 3D muon propagation code that takes into account energy losses due to bremsstrahlung, pair production, 
inelastic scattering and ionization. It simulates also the angular and lateral deflections due to multiple scatter- 
ing. In order to sample the showers on the detector surface, a cylindrical volume, surrounding the sensitive 
volume of the detector, is defined around the PMT array, the can. Its dimensions are obtained adding 2.5 x Xabs 
to the maximum distance of PMTs from the centre of gravity of the PMT array, where Xabs — 55 to is the 
water absorption length at the ANTARES site. The dimensions of the can used in this work are : radius ~ 240 
m; height ~ 590 m. 

For each shower a random point is generated on the surface of a cylinder that is larger than the can, in order to 
take into account the lateral extension of the shower Starting from this impact point, all the possible intersec- 
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tion points of each muon with the can surface are evaluated. Each particle is then propagated along the path 
between the sea level and the detector using the MUSIC code, requiring a minimum residual energy of 20 GeV 
at the surface of the can. 

Figure 1 shows the energy spectrum and the angular distribution, normalised to one day of live time, of the 
muons reaching the total can surface. The continuous line, corresponding to the HEMAS simulation, and the 
dotted line, relative to the CORSIKA simulation, are almost perfectly coincident. The total number of muons 
is about 2.7 x 10'' with a difference smaller than 4% between the two productions. 

The third step of the simulation concerns the detector response and is performed using a Geant3 j9| based 
Monte Carlo code. An optical background rate of 60 kHz per PMT has been added. This is the rate expected 
from the presence of A' in sea water and has been measured in situ LlOi . 



3. Track reconstruction. 

Finally, the hits registered on the PMTs have been processed in order to reconstruct the track direction. Several 
strategies have been developed in ANTARES for track reconstruction. In this work two approaches have 
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Figure 2. Angular distribution of reconstructed muon tracks with the 1111 strategy for: HEMAS sample (left); CORSIKA 
sample (right). The continuous lines refer to the reconstructed angles, the dotted lines to the original direction of events. 



been considered, both of them based on the calculation of a likelihood function, starting from the recorded 
information of the arrival time and of the amplitude of hits on each PMT. Its maximum is used to obtain an 
estimate of the muon track parameters. 

The first strategy is described in detail in II II . The results obtained when applied to the muon atmospheric 
samples from the HEMAS and CORSIKA simulations are shown in figure 2. The continuous lines refer to 
the reconstructed angles, the dotted lines to the original direction of events. About 2.6 x lO'' downward going 
muon tracks per day are reconstructed. Also some fake upward going events are reconstructed: 8 ± 2 per day 
for HEMAS and 10 ± 2 for CORSIKA. About 40% of them correspond to muons arriving in a bundle. 
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The other strategy used for reconstruction is described in 1 12|. It is optimized to reconstruct upward going 
tracks and its efficiency in reconstructing downward going atmospheric muons is poor with respect to the first 
strategy. On the other hand, it has a better capabiHty to reject misreconstructed upward going muons (no 
upward going tracks have been reconstructed, when appHed both to the HEMAS and the CORSIKA samples). 

4. Conclusions 

A full simulation of the atmospheric downward-going muon flux has been performed using two atmospheric 
shower simulation programs, HEMAS and CORSIKA. No significant differences have been found at the de- 
tector level (see fig. 1) between the two samples. 

An evaluation of the contamination due to the atmospheric downward-going muons erroneously reconstructed 
as upward going is presented, using two reconstruction strategies available in the ANTARES Collaboration. 
One of them 1 1 1 1 shows a better efficiency in downward track reconstruction, while the other one 1 12 1 rejects 
all upward going reconstructed tracks. The choice of the most appropriate strategy for reconstruction strongly 
depends on the characteristics of the analysis to be performed (diffuse neutrinos, pointlike neutrino sources, 
atmospheric muons etc.). An excellent agreement between the two samples is obtained also after the recon- 
struction step. 

The results are very encouraging. Work is in progress to improve their statistical significance, increasing the 
number of simulated events. 
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